Extracting the hydrodynamic resistance of droplets from their behavior in 

microchannel networks 
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The overall traffic of droplets in a network of microfluidic channels is strongly influenced by the 
liquid properties of the moving droplets. In particular, the effective hydrodynamic resistance of 
individual droplets plays a key role in their global behavior. We here propose two simple and low- 
cost experimental methods for measuring this parameter by analyzing the dynamics of a regular 
sequence of droplets injected into an "asymmetric loop" network. The choice of a droplet taking 
either route through the loop is influenced by the presence of previous droplets which modulate the 
hydrodynamic resistance of the branches they are sitting in. We propose to extract the effective 
resistance of a droplet from easily observable time series, namely from the choices the droplets make 
at junctions and from the inter-droplet distances. This becomes possible when utilizing a recently 
proposed theoretical model, based on a number of simplifying assumptions. We here present several 
sets of measurements of the hydrodynamic resistance of droplets, expressed in terms of a "resistance 
length". The aim is twofold, (1) to reveal its dependence on a number of parameters, such as the 
viscosity, the volume of droplets, their velocity as well as the spacing between them. At the same 
time (2), by using a standard measurement technique, we compare the limitations of the proposed 
methods. As an important result of this comparison we obtain the range of validity of the simplifying 
assumptions made in the theoretical model. 



I. INTRODUCTION 

Droplets in microfluidic channels are widely used as 
microreactors for chemical analysis, for chemical kinetics 
studies, or, for examplejas a tool for continuous nanopar- 
ticle synthesis [H, i, M 11 H- Propelled by the necessity 
to accommodate parallel processes and to sample or sort 
droplet populations, the demands on microfluidic net- 
works to implement certain procedural targets keep on 
increasing. One strategy towards the parallel processing 
of many droplets are self-regulated channel networks in 
which it is not necessary to control the route of individual 
droplets by external means (valves) but in which a well- 
chosen channel design leads to a autonomous behavior of 
the droplet train [2,0,0, The advantage of this ap- 
proach is evident as the channel geometries stay simple, 
and less external devices are required. A drawback of the 
approach is an increased complexity in the understanding 
of the droplet dynamics, making it more difflcult to flnd 
robust implementations of given tasks [1, [l^, lU • Gen- 
eral features of self-regulated droplet dynamics are the 
global coupling of the dynamics, meaning that nearly all 
droplets influence each other, and the strong nonlinear- 
ity in the equations describing the droplet positions as a 
function of time. 

The non-local and non-linear droplet dynamics has 
been described in Ref. [ll| by a simplified theoretical 
model. In this model, each droplet, sitting in one con- 
necting channel, increases the hydrodynamic resistance 
of this channel by a given constant amount. By its pres- 
ence it changes the flowrate in the channel, and in turn it 
changes also the flowrates in possibly all other channels of 
the network. When a droplet arrives at a junction, a rule 
is required to decide which of the possible exits is to be 
taken. It will here be assumed that it takes the channel 



with the highest instantaneous flowrate. The simplified 
model of Ref. ll| was able to address several intriguing 
properties of different network layouts, such as the re- 
versibility of the droplet dynamics and the influence of 
network symmetry. The model was proposed as a simple 
tool for flnding robust dynamical behavior and to quan- 
tify its response to changes in the driving parameters and 
the geometrical parameters of the network. While apply- 
ing the model calculations to existing devices, it turned 
out that the most important parameter for robust behav- 
ior is the effective hydrodynamic resistance of droplets. 
Compared to its important role in the overall dynamics 
of droplets in a network, this property has not attracted 
sufRcient interest, at least for droplets of a viscous fluid. 
As a first guess, it was therefore assumed to be a constant 
in the model. 

In this paper, we present results from three methods for 
measuring the hydrodynamic resistance of droplets being 
transported in microchannels: A standard method that 
employs a pressure sensor as a reference, and two new 
self-regulated low-cost methods which do not require la- 
borious external equipment. These novel methods follow 
the idea of the abovementioned self-regulated networks. 
We here propose to make use of the overall droplet dy- 
namics, which is easily observable with a video device, 
in order to extract effective properties such as the hydro- 
dynamic resistance. As we make use of self-regulating 
droplet dynamics, which requires a theoretical model for 
its interpretation, the target here is twofold: As long as 
the assumptions of the model from Ref. [11] hold, the 
property extracted from the measurement is really the hy- 
drodynamic resistance of the droplets. At the same time 
we verify these very assumptions by checking the consis- 
tency of the outcome and by directly measuring the hydro- 
dynamic resistance of droplets with the standard pressure 
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FIG. 1: Schematic pictures of the devices: (a) A straight 
channel with two pressure taps for the method using the pres- 
sure sensor, and (b) an asymmetric circular loop. Droplets 
are generated where the first oil and water/glycerol inlets 
meet. The third inlet allows to adjust the distance between 
the droplets. Channels and inlets have a rectangular cross- 
section, 500 /im wide and 300 fim high. Only the two channels 
to the pressure sensors are narrower (100 /im wide). 
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FIG. 2: Time series of the pressure difl^erence along a straight 
channel, as measured with the pressure sensor. The droplet 
input has been changed subsequently: First (1) without 
droplets, then (2) with them, and (3) last without droplets. 
The total flow rate has been kept constant. The dotted lines 
are guides for the eye. Flow rates of the oil and the wa- 
ter inlets in the droplet generator are respectively 8/il/min 
and 7 fjl/ min. The second oil inlet allows to adjust the dis- 
tance and the velocity of droplets, here 70/il/min. The pa- 
rameters of the droplets are rjd = 0.001 Pas, Qd ~ 85 nl, 
Vtot = 7.5mm/s, A = 6 mm. 



sensor. 

The paper is organized as follows: First, in Sec. [Til we 
explain the principles of the methods the specific advan- 
tages and limitations of each of them: First, the well es- 
tablished method based on the use of an external pressure 
sensor, and then, the principle of the new methods rely- 
ing on the flow of a regular train of droplets in a "loop" 
device. We distinguish between two different sizes of the 
loop, giving different weights to the underlying assump- 
tions on which the interpretation is based. In Sec. Illll we 
describe the experimental procedure we used for creating 
the device. Section llVl contains the experimental findings 
of the two variants of the loop method and of the method 
using a pressure sensor. There, we discuss in detail the 
implication of the measurements, being in favor of the 
simple model or putting its assumptions into question. 



II. PRINCIPLES OF THE METHODS 
A. The pressure sensor method 

Experiments were conducted using a differential pres- 
sure sensor as described by by Adzima et al. [T5|. At a 
distance of 127 mm, two pressure taps were connected to 
a long straight channel, as depicted in Fig. [1^. At each 
of them the pressure is measured by a pressure sensor 
(Validyne P55D), equipped with a 1.25 psi-membrane. 
By comparing the pressure differences in two situations, 
namely with and without droplets — while keeping the to- 
tal flowrate Q the same — we infer the additional pressure 
caused by the presence of N droplets between the taps. 
The hydrodynamic resistance Rd of each droplet is then 



given by 



Rd 



AP* - AP 



(1) 



where AP* and AP are the pressure differences between 
the taps, with and without droplets respectively [l2|. In 
the discussion below it will prove convenient to express 
hydrodynamic resistances in terms of lengths. This is 
natural for channels without droplets where the hydro- 
dynamic resistance Ri of a channel i is known to scale 
linearly with its length L^, 
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Here, w and h are the channel width and height, respec- 
tively; 77$ is the viscosity of the continuous phase; and 
the dimensionless parameter A depends on the shape of 
the channel cross-section, which is rectangular here p^ . 
We apply the same factor of linearity from Eq. ((J) also 
to the effective resistances Rd of the droplets in order to 
define their resistance length Ld^ 
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The pressure sensor is used below in Sec.|IV]to compare 
the results of the alternative passive methods presented 
in the following subsections. It is important to note that 
this equipment is quite costly (several hundreds of dol- 
lars) and that it presents a number of limitations in terms 



FIG. 3: Snapshots of the small-loop device, fed by a regular 
sequence of droplets: (a) all droplets take the short route, 
the incoming distance is large, A > A,nax- (b) alternating 
droplet decisions due to interactions between droplets, A < 
Amax- The arrows indicate the flow direction. Parameters are 
Ls = 4.25 mm, Le/Ls = 1.18. 



of response time and ease of use: The smaller the probed 
channel, the larger is the equilibration time. In the case 
of our channels, which are quite large (300 /im high and 
500 /im wide) , this equilibration time was found to be of 
the order of ten minutes. Figure [5] shows the pressure 
difference at the two taps as a function of time, exhibit- 
ing the slow relaxation of the signal at the beginning and 
at the ending of a droplet series. The following passive 
methods allow to get rid of these constraints. 



B. The loop devices 

We here propose two variants of a new method which 
allows to extract the hydrodynamic resistance of droplets 
from the observed behavior of regular droplet trains. 
This approach has the advantage over the standard pres- 
sure sensor that it is cheap and does not require taps 
which pose an additional source of perturbation of the 
flow. Furthermore, it can be applied without problems to 
very small channels, where the pressure sensor technique 
becomes inapplicable. Of course, our method has its lim- 
itations: The main limitations come from the hypotheses 
made in the theoretical model on which the interpreta- 
tion of recorded droplet behavior is based. It will turn 
out in Sec. IIVI that these assumptions are not always 
satisfied. Nevertheless, in those regions where they do 
hold, the new method presents a real alternative to the 
pressure sensor. 

We use a very simple network of microfluidic channels 
as depicted in Fig. (TJa. Disregarding the part in which 
droplets are generated, it consists of two long channels 
with an asymmetric circular loop in-between. A periodic 
sequence of droplets is generated and injected into this 
loop from one side. The droplets then take either the 
short or the long branch of the loop, depending on how 
many droplets the respective branches contain: For the 
first droplet, the shorter branch is certainly favorable. 
It has a higher flowrate and a smaller hydrodynamic re- 
sistance. The arriving first droplet then increases the 
hydrodynamic resistance of this shorter branch, possi- 
bly so much that the following droplet takes the longer 
route. Upon continuation, a complicated sequence of 
droplet "decisions" is generated, exhibiting complex pe- 
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FIG. 4: Snapshots of the large-loop device: (a) the arriving 
droplet chooses the long branch, corresponding to the rela- 
tion (|14[) . here Re + 3Rd < Rs + 5Rd- (b) the next arriving 
droplet takes the short route, Ra + iRd < Re + SRd- The 
dotted arrows indicate the direction chosen by the droplets. 
Parameters are Ls — 37 mm, Le/Ls — 1.16. 

riodicity. The same network topology has been used 
in Refs. [l^, |Tl| to investigate the reversibility of the 
droplet dynamics. The droplet traffic in the network was 
recorded by a charge-coupled device (CCD) camera and 
processed with a homemade Matlab program. Such, we 
extracted from the movies the velocities of all droplets, 
the number of droplets in each channel, and the tempo- 
ral and spatial periods of the incoming droplet sequence. 
These properties were then plugged into the theoretical 
model in order to extract the hydrodynamic resistances, 
according to the two methods described in the following: 
The flow through the devices is governed by flow equa- 
tions which resemble Kirchhoff's law for electric circuits. 
At each junction, the sum of incoming and outgoing 
flowrates must vanish. This means for the loop device 

Qtot ^ Qs + Qi, (5) 

with Qtot the total flowrate in the input/output chan- 
nels, and with Qs and Qi the flowrates in the short and 
long branches of the loop. Here and in the following, the 
indices s and £ stand for the short and the long chan- 
nel, respectively. The second set of equations provides 
the relationship between the pressure difference along a 
channel and the flowrate passing through it, 

{AP), = R,Q,. (6) 
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FIG. 5: Experimental test for the linear relation between 
droplet velocity Vi and the flowrate Qi (Assumption A). 
The data were acquired in a straight channel. The param- 
eters are r/d = 0.76 Pas, 77$ — 0.02 Pas, A = 1.7 to 2.8mm, 
D,d = 55 nl. 



At the moment, this equation is nothing but the defini- 
tion of the total hydrodynamic resistance Ri of channel i. 
Without any droplets, it is known to provide a linear re- 
lation, in that the resistance does not depend on the total 
flowrate. 

Before describing the dynamics of droplets in the two 
geometries we used, let us summarize the assumptions 
which led to the simplified model and which will be in- 
vestigated experimentally in the following sections: 

(A) All droplets in a channel i have the same velocity 
velocities Vi. Its value is proportional to the total 
flowrate Qi in the channel (oil -I- water/glycerol), 
with a constant factor of proportionality. 



Vi = — 



(7) 



Here, S is the cross-section which does not change. 
(3 was found to be approximately 1.6 in our exper- 
iments. 

(B) Each droplet adds the same amount Rd to the total 
resistance Ri of a channel i, such that 



Ri — Ri + TiiRd, 



(8) 



with rii the number of droplets in the channel, and 
with Ri its resistance in the absence of droplets. 

(C) The resistance Rd of droplets should be a con- 
stant. In particular, Rd should not depend on the 
flowrate Qi (nor on Vi). 

(D) At junctions, droplets take the route with highest 
flowrate leading away from the junction. 

Some comments on the assumptions and their expected 
validity are in order at this point: The linearity of ve- 
locity and total flowrate (A) has been observed under 



experimental conditions, both in Fig. [5] and in Ref. p^ . 
Assumption (B) will prove false if there is cooperative 
behavior between the droplets. Such effects will not oc- 
cur as long as the droplets are sufficiently distant. We 
will quantify the necessary distance in the experimental 
section. Assumption (B) will also fail if the droplets have 
different sizes or if they are smaller than the channel and 
are thus not forced to be centered. Assumption (C) is 
nothing else than the linearity of Eq. ([6]). In fact, we 
will present deviations which are important for small ve- 
locities, but not for larger ones. Also in case of long 
plug- like droplets, a nonlinearity in Eq. ([H]), and there- 
fore a velocity-dependence of Rd must be expected. Such 
effects have been reported for bubbles [H, [H, [13, Ell . By 
construction of our device, that is by virtue of its sim- 
ple geometry, we do not require assumption (C) for the 
method to be valid. Finally, assumption (D) will not be 
questioned here. The precise criterion, whether it is the 
highest flowrate, the highest pressure drop, or the largest 
velocity which counts, does not seem to be important in 
our case as we use locally symmetric junctions and the 
same channel cross-sections everywhere. 



1. The small-loop device 

The first variant of the proposed device comprises a 
small loop, as shown in Fig. [31 The two branches of 
the loop are approximately ten times longer than they 
are wide. The idea of the method is to send in a pe- 
riodic sequence of droplets and to vary slowly their ini- 
tial distance. If one starts with quite a large distance, 
the droplets will all take the short route and leave the 
loop before the next one arrives, see Fig. [3^. For a con- 
tinuously reduced initial distance A, one finds a sudden 
change in the droplet behavior at a specific distance Amax- 
The droplets start to influence each other: The following 
droplet will take the long route because the previous one 
still enhances the resistance of the short branch. This sit- 
uation is depicted in Fig.[5lD. Of course, in order to make 
the method work, the (yet unknown) hydrodynamic re- 
sistance of a single droplet must be larger than the differ- 
ence in resistances of the channels when no droplets pop- 
ulate them. The applicability can therefore be adjusted 
by choosing the length ratio LijLs- The hydrodynamic 
resistance of the droplet is determined from the value 
of the initial distance Amax, at which the droplets stop 
taking the entire short branch, by applying the simple 
Eq. (jl3p below, which we will derive now. 

The initial distance A,nax is obtained as the distance 
which an approaching droplet can travel while the previ- 
ous one traverses the short branch. The time T which a 
droplet stays in the short branch is expressed as follows, 
when assumption (A) is used. 



T ■ 



PQs 



(9) 



where L^, Vs, and Qs are the length, the velocity and 
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the flow rate in the short arm. In the same time span, 
another droplet in the incoming channel covers the length 



(10) 



which leads to the elimination of the parameters /3 and 5, 
and which leaves us with 



^max 



Qs 



= 1 



(11) 



Now, we make use of assumption (B), together with 
Eq. As the pressure difference is the same in both 
branches of the loop, we obtain for the ratio of their 
flowrates 



Qe 



Ri 



Le 



(12) 



Here, it proves convenient to express the hydrodynamic 
resistances in terms of lengths, as defined in Eq. 
Combining the last two equations, we find the result al- 
ready given in Ref. 11], 



An 
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Ld 
L/ 



(13) 



which can easily be resolved for Ld- Note that we did 
not make use of assumption (C) in this short derivation. 
Equation therefore holds for any non-linear relation- 
ship between Qs and the resistance Rd{Qs) of the droplet. 
It is the simple geometry of the device which allows to 
circumvent the simplifying assumption (C), which is not 
possible in extended networks. Note also the nice prop- 
erty of this formula that the hydrodynamic resistance 
length of the droplet is only given by geometrical prop- 
erties, namely the two lengths of the branches, and by 
the experimentally determined initial distance Amax be- 
tween droplets. This means that once the device has been 
elaborated, the measure of the droplet resistance is sim- 
ply obtained from observing the droplets taking different 
routes through the loop. 

Note also that the determination of the droplet resis- 
tance from this device relies strongly on the assump- 
tion (A), and only weakly on the assumption (D). As- 
sumption (C) is not used at all, and assumption (B) 
is unimportant, since we treat one droplet per channel 
(n = 1) which is the trivial case for (B). 

The limitation of the here described small-loop method 
is that the two branches are quite short. A violation 
of the above derivation might occur, in that the resis- 
tance Rd depends on the instantaneous distances between 
droplet and the two junctions, both of which perturb the 
flow in a non-trivial manner, leading to an effective hy- 
drodynamic interaction between them. In this aspect, 
the smaller the channel cross-section is, the better the 
small-loop device works. This will lead to a more uni- 
form flow and at the same time increase the resistance of 
droplets. 



2. The large loop deviee 

We propose a second variant of the loop device which 
allows to determine the hydrodynamic resistance of 
droplets with a different weighting of the assumptions. 
This time, the loop is much larger compared to the pre- 
vious geometry, see Fig. 01 Like for the small loop device, 
a regular train of droplets is injected into the main chan- 
nel. This time, there can be several to many droplets 
in each branch of the loop, and their resistance is de- 
termined by counting the droplets in each branch at ev- 
ery time a new droplet has made its decision. Thus, the 
droplet dynamics depends on the model assumptions in a 
different way: According to assumption (D) each droplet 
entering the loop systematically goes into the branch of 
lower hydrodynamic resistance. And, since each droplet 
adds a certain resistance to the channel (assumption B), 
it tends to equilibrate the hydrodynamic resistances of 
the two branches of the loop. The same principle has 
been previously used in the analysis of the hydrodynamic 
resistance of bubbles ^1^, however only in the case of a 
symmetric loop and without the complementary device 
with a small loop. 

During the run of the experiment, the total hydrody- 
namic resistances of the two branches change as functions 
of the numbers of droplets they already contain. We thus 
have to track the relation between these two resistances. 



Rs + UsRd ^ Ri + niRd, 



(14) 



where both comparison operators are found. Each time 
a droplet has chosen a route, we gain some information 
about which branch of the loop has had the smaller to- 
tal resistance. In fact, each time a droplet chooses the 
longer branch, we obtain a lower bound for Rd] the 
short branch provides an upper bound. When resolving 
Eq. (jl4p for Rd, rewritten in terms of the corresponding 
length Ld, we obtain these bounds as 



Lt — Ls 
Hs - ni 



< Ld < 



^long 



Le — Ls 
Hs - ne 



(15) 



short 



The subscripts ^^—^long" and ^^^shorf^ indicate that the 
specified bound applies only if the long/short branch 
has been chosen. Observing many droplets will succes- 
sively narrow the bounds on Ld, leaving finally an interval 
which depends on both the initial droplet spacing A and 
on the length ratio Li/Lg- A large loop which is nearly 
symmetric leads to a good accuracy of the method. Fig- 
ure m shows two examples of how Eq. (fT5)) is applied, 
leading in Fig. 2^ to a lower bound, and in Fig. 3)3 to an 
upper bound. 

In terms of the assumptions required for the analysis, 
this large-loop method appears to be complementary to 
the small-loop variant. We here make heavy use of as- 
sumption (B) in that we assume the droplets to add up 
their individual resistances. On the contrary, the small- 
loop proposition was mainly based on assumption (A) 
which is not really required here, because the time the 
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droplets stay in in both branches of the loop are not used. 
Also the role of the droplet velocity is different. While the 
flowrates in the two branches of the small-loop device can 
be very different, they are automatically balanced here — 
at least their average values, with fluctuations due to the 
fluctuating numbers of droplets. 

III. EXPERIMENTAL DETAILS 

We made the microfluidic devices out of poly(dimethyl- 
siloxane) (PDMS) using standard pho to-lithographic and 
soft-lithographic techniques [13, l2l| . A transparency 
mask of 20 000 dpi resolution (CAD/Art Services Inc., 
California) was used in 1:1 contact photolithography with 
SU-8 photoresist (MicroChem, Newton, MA) to generate 
a negative "master" consisting of patterned SU-8 pho- 
toresist on a silicon wafer with a flat surface. Positive 
replicas were formed by molding a 5:1 mixture of PDMS 
with a temperature-active cross-linker against the mas- 
ter. After curing the PDMS layer for 20 minutes at 65° C, 
it was peeled off the silicon wafer, and inlet and outlet 
holes were punched. A glass slide, serving as the rigid 
substrate of the device, was covered by a thin layer of 
a 20:1 mixture of PDMS/cross- linker and cured during 
40 minutes at 65° C. The PDMS surfaces were then ac- 
tivated with an UVO-cleaner (Jelight 144AX) [111 for 
40 s, and both PDMS surfaces were immediately brought 
together. An irreversible seal was formed between the 
PDMS surface and the glass substrate. 

The resulting rectangular microchannels are 300 /zm 
high and 500 fim wide. Droplets were fabricated from 
a droplet generator with a cross-junction (Fig. [T|): In- 
lets for different liquids, one of oil (viscosity 77$ = 20 cP) 
and one for a glycerol/water mixture, lead into a channel 
located downstream the junction. A colorant (Bromoth- 
ymol blue) was added to this aqueous phase in order to 
increase the contrast of the droplet. External connec- 
tions to syringes (Hamilton, 2.5-10 ml) were made us- 
ing polyethylene (PE20) medical tubing (outer diameter: 
1.09 mm). The syringes were contracted by three pro- 
grammable syringe pumps (Cetoni, Nemesys) controlled 
by a computer. An additional oil inlet further down- 
stream served to vary the spatial period between droplets 
without modifying the volumes of the droplets. Their 
volumes were inferred from all three incoming flow rates 
and the measured droplet distance. 



IV. RESULTS AND DISCUSSION 

In this section we present several sets of measurements 
in which different parameters are varied. We measure 
their influence on the hydrodynamic resistance of the 
droplets, in order to see to what extent the assumptions 
(A) to (D) may be applied. The varied parameters are 
the volume of the droplets, their velocity, their vis- 
cosity rjd and their initial distance A. All these parame- 



ters influence the hydrodynamic resistance of droplets — 
which, strictly speaking, invalidates assumptions (B) 
and (C). Nevertheless, we will show that reliable results 
can be obtained with the newly proposed methods. Their 
reliability is enforced by the comparison with the results 
from the pressure sensor as a reference. The first two 
presented results coincide well and show that the pro- 
posed small-loop and large-loop devices can be used to 
determine the hydrodynamic resistance of droplets. 

From the remaining results we learn that the assump- 
tions underlying these methods are not fulfilled for all 
driving parameters. On first sight, this appears as bad 
news for the methods we propose. However, we also find 
that there is always a parameter regime in which the as- 
sumptions are indeed valid. Concretely, we find that both 
the droplet velocity and the inter-droplet distance must 
be kept sufficiently large. If we ensure that the velocity 
does not drop below a critical value, while varying for 
example the distance, the result is perfectly valid. We 
may use such a situation to extract the hydrodynamic 
resistance of the droplets. The results given in this sec- 
tion therefore present both a direct verification of the 
assumptions of the employed model and a measurement 
of the hydrodynamic resistance. 

A. Validity of the assumption (A) 

Independently of the parameter variations. Fig. [S] 
shows that the assumption (A) holds in the performed 
experiments. 



B. Influence of the droplet volume 

The influence of the volume of droplets on their hydro- 
dynamic resistance has been investigated with the pres- 
sure sensor and with the small-loop device. The results 
in Fig. [6] depict a strong increase of the resistance with 
the volume, which is confirmed by the measurements of 
the pressure sensor. Both methods yield coinciding val- 
ues for the resistance. The increase can be understood 
by the fact that with growing volume the droplets fill the 
channel more closely and become longer. Similar effects 
have been reported in Refs. [H, [l3|. Droplets of 48 nl 
are nearly spherical when looked at from above. Larger 
droplets are more elongated, giving rise to thin films of 
oil between the droplets and the channel walls. 

The viscous dissipation in these thin films is one of 
three possible reason for the enhancement of hydrody- 
namic resistance due to the presence of a droplet: Gener- 
ally speaking, each droplet perturbs the surrounding fiow 
field, which would otherwise be an axial Poiscuille-type 
flow, see Fig. [7] and Ref. [2^. The kinematic bound- 
ary condition at the boundary of the droplet enforces 
the presence of several counter-rotating eddies inside and 
outside the droplet. Figure [7] contains only those eddies 
which are necessary to ensure a continuous velocity field. 
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FIG. 6: The resistance length Ld as a function of the droplet 
volume: measured with the small-loop device and with the 
pressure sensor. The parameters are rjd = 1.4 Pas, 77$ — 
0.02 Pas, La = 4.25mm, Le/La = 1.18, vtot ~ 6mm/s, and 
A > 9 mm. 
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FIG. 7: The topology of the counter-rotating eddy pairs in 
front, behind and within the droplets moving in a channel. 
The velocity streamlines are those seen from a frame co- 
moving with the droplets: The channel walls move to the left. 
Depending on the viscosity ratio, between one and three eddy 
pairs can be found in the droplets. More eddies are possible, 
depending on the shape of the droplets, thus on the surface 
tension. 



FIG. 8: The resistance length Ld (large-loop method and 
pressure sensor) as a function of the droplet velocity Wtot. 
Inset: length of the droplet (cap to cap). The parameters 
are rjd ~ O.OlPas, = 0.02 Pas, La = 53mm, Lt/La — 1.5, 
Q,d ~ 95 nl. The distance varies as A from 5 mm to 15 mm 
from small to large velocities. 

distance between droplets. 

In the experiment of Fig. [H] the range of droplet vol- 
umes is limited, as large droplets tend to break or deform 
and wait in the incoming junction of the loop: The wait- 
ing time in the junction causes an overestimation of Ld by 
Eq. ini). Note that in the small- loop model of Sec. Ill B II 
the new route is chosen instantaneously and does not 
take any waiting time into account. There is also a lower 
bound for the applicable droplet volume, given by the 
channel cross-section. As soon as the droplets have no 
contact with the channel walls, their resistance strongly 
decreases. However, we require their resistance to be 
large enough to overgrow the asymmetry of the loop. 

C. Influence of the total flowrate 



In particular, the tangent velocity has to be continuous 
at the droplet boundary. Due to their different viscosi- 
ties, the eddies inside and outside the droplet may con- 
tribute differently to the total hydrodynamic resistance of 
a droplet. A third contribution comes from the thin fluid 
film between droplet interface and channel wall. Its re- 
duced thickness leads to large velocity gradients and thus 
also to viscous dissipation. As the droplet resistances at 
different volumes are compared, it is not only the thin- 
film contribution which changes. Together with the more 
elongated shape of the interfaces, also the precise shape of 
the interior eddies change. The viscous dissipation within 
the droplet strongly depends on the details of the eddies, 
in particular on the ratio between the round backflows 
near the caps of the droplet and the rather "Poiseuille- 
like" parts in the interior of a long droplet. At this point, 
nothing can be said about the relative importance of the 
viscous dissipation in the thin film and in the eddies. We 
will return to this point below, in the discussion of the 



According to assumption (C) the hydrodynamic re- 
sistance length Ld should neither depend on the total 
flowrate in the channel, nor on the velocity of the droplet. 
In fact, we found that this is the case only for velocities 
higher than approximately Vtot — 5 mm/s, as Fig.[8]indi- 
cates. At small droplet velocities we find very large hy- 
drodynamic resistances. The droplet resistance decreases 
with increasing velocity, until only a weak dependence is 
observed beyond Vtot = 5 mm/s. The data plotted in 
Fig. [H] have been obtained from the large-loop device and 
are in good agreement with the reference measurement. 
It is important to stress this agreement, as it shows that 
the large-loop device indeed measures the droplet resis- 
tances, and that it does not depend on assumption (C), 
as was already explained at the end of Sec. IIIB 21 

The velocity value vtot in the measurement is the one 
in the incoming channel, which is systematically differ- 
ent from the true velocity at which the droplets move. In 
the large-loop device, however, we find automatically bal- 
anced flowrates and therefore also balanced velocities in 
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FIG. 9: The resistance length Ld (smaU-loop method) as a 
function of the droplet viscosity. The parameters are rj^ — 
0.02 Pas, Ls = 4.35mm, Le/Ls = 1.18, vtot = 9mm/s, A > 
9 mm, and Qa = 65 nl. 



the two branches, which allows to say that both are half 
the incoming velocity, Vg ^ Vi ^ utot/2, at least on aver- 
age. That we obtain the true velocity immediately from 
a measurement in the straight part of the channel is the 
reason why we performed this experiment with the large- 
loop device and not with the small-loop device, where the 
velocity ratios depend on the droplet resistances, unless 
the droplet velocity is sufficiently high. 

The measurements in Fig.[5]have been performed using 
droplets of the same volume and the same composition. 
The velocity was varied by injecting more oil between the 
droplets. Together with the total flowrate, thus also the 
initial spacing between the droplets was varied between 
5 mm and 15 mm. The origin of the strong increase of re- 
sistance at small velocities is unclear at the moment. We 
found a correlation with the shape of the droplets (inset 
of Fig. [5]), indicating a systematic lengthening of the fast 
flowing droplets. This increases the thickness of the thin 
oil film between the droplet and the walls and reduces the 
total resistance of the droplets ■ A dynamical theory 
which would allow to fit the droplet resistances given in 
Fig. [H] is not known to us. In fact, we tried to fit the data 
to existing formulae for the resistance as a function of the 
droplet size and of the velocity: In Ref. [l^ a proposition 
containing a scaling with u^/^ is put forward for bubbles. 
In Ref. ^] the scaling is rather v^^"^. Neither formula 
yielded a remotely satisfactorily result for the data in 
Fig. [51 The details of the surrounding flow and of the in- 
ternal circulation in the droplet are evidently much more 
important than in the case of bubbles, where a simple 
scaling approach seems to work [l^, [l^, [2^ . 



D. Influence of the droplet viscosity 

Another set of experiments with the small-loop de- 
vice explores the influence of the viscosity of the glyc- 
erol/water mixture inside the droplets on their resistance. 



The viscosity was changed by varying the concentrations 
of glycerol and of water in the mixture, the continuous 
phase and the velocity were kept constant. Figure [S] de- 
picts the resulting resistance length Ld as a function of 
the fluid viscosity in the droplet. The values of remain 
approximately constant as long as the droplet viscosity 
is lower than the viscosity of the surrounding phase. Be- 
yond this value, they exhibit a strong increase which can 
be attributed to the enhanced viscous dissipation by the 
internal circulation in the droplets. Similar effects have 
been observed in Ref. [13] ■ The agreement between the 
results of the small-loop device and the pressure sensor 
is mediocre. We do not have a convincing explanation 
for the deviations at the moment. Nevertheless, the two 
measurements exhibit the same ascending trend. As de- 
scribed above, the velocity of the droplets were quite fast, 
in order to stay within the part of Fig. [8] where the ve- 
locity does not strongly influence the resistance length. 
We chose this regime in order to clearly separate the re- 
spective influence of the viscosity from others. 

The qualitative form of the flow fields are depicted in 
the sketch of the counter-rotating eddies in Fig. [71 The 
two panels of the figure present two different topologies 
for the two extreme cases of the ratio rjd/r],^. These two 
topologies are expected because the fluid with higher vis- 
cosity should avoid to produce several eddies, each of 
which increases the viscous dissipation. Despite the fact 
that there are two clearly distinguished regimes of the 
function Ld{rid) in Fig. [51 in one of which it is nearly con- 
stant, in the other fast growing, we flnd a rather smooth 
change from one regime to the other, and not an abrupt 
one at the point of equal viscosity rjd = rj^. The curve 
may further be extrapolated to an infinite droplet viscos- 
ity, in which the result of a nearly filling rigid object in 
a channel is reached. The resistance of such an object 
depends solely on its filling property, that is on the thin 
film between object and walls. It is then only the ratio of 
surface tension and exterior viscosity (exterior Capillary 
number) which determines the thickness of the thin film. 



E. Influence of the initial inter-droplet distances 

We now discuss possible cooperative effects between 
adjacent droplets and thus explore the validity of as- 
sumption (B). Of course, this can only be undertaken 
by the large-loop device, since the small loop contains 
not more than one droplet. Upon varying the initial dis- 
tance between droplets, we find indeed a critical droplet 
spacing below which assumption (B) does not hold any 
more. Figure [lOl visualizes this result, as measured by 
the large-loop device. For small distances, an increase 
of the droplet resistance as a function of the distance 
can be seen. Beyond approximately 6 mm we find an ap- 
parent plateau. The large uncertainties in this part of 
the data are caused by the small number of droplets in 
the branches of the loop. For any given maximal size of 
the large loop results a maximal reasonable droplet dis- 
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FIG. 10: The resistance length Ld (large-loop method) as a 
function of the initial inter-droplet distance A. Inset: vari- 
ation of the droplet velocity vtot with A. The parameters 
are rjd = 1.4Pas, ?7$ ~ 0.02 Pas, Ls = 36mm, Li/Ls — 1.6, 
fid = 65 nl. 



tance. The uncertainties can, however, be reduced either 
by using a ratio L^/La closer to unity, or by investigat- 
ing narrower channels. In the latter case, the droplet 
distance up to which cooperative effects occur, should be 
decreased together with the width/height of the channel. 
The necessity to work at distances larger than 6 mm in 
order to guarantee the validity of assumption (B), has al- 
ready been taken into account in the previously described 
measurements. 

The decrease of resistance at distances smaller than 
6 mm gives evidence to a collective motion of the droplet. 
This means that a pair or series of droplets can be dis- 
placed more easily than the same number of individual 
droplets. This observation sheds light on the relative 
importance of the three different contributions to the hy- 
drodynamic resistance, as mentioned above in the discus- 
sion of Fig. [71 In the literature on (inviscid) bubbles, the 
contribution from the thin film has been considered as 
the primary reason for hydrodynamic resistance [isl [T^ . 
The decrease of resistance observed in the measurements 
of Fig. [TO] proves that in our case the thin-film contribu- 
tion cannot be dominant but that the eddies outside of 
the droplet contribute equally. From Fig. [TOj we can es- 
timate the size of the eddies. The incoming distances of 
6 mm correspond to a minimal distance between droplets 
of around 3 mm in the branches where the droplet se- 
quences are not regular anymore. The diameter of the 
exterior eddies are then around 1 mm, thus two times 
the width of the channel. 

Note that for experimental reasons, the three rele- 
vant parameters, namely droplet velocity, droplet spacing 
and the droplet volume cannot all be independently con- 
trolled. Therefore, this set of experiments has not been 
performed at constant flowrate (inset on Fig. [TO)) . There 
is thus a hidden influence of the droplet velocity on Ld, 
which should pronounce the dependence on A even more. 

The full function Ld{X,vtot,^d,'nd,V'S>) can only be 
fully displayed in a high-dimensional plot. A non-trivial 
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FIG. 11: Combined influence of two parameters: The re- 
sistance length Ld (all three methods) as a function of the 
droplet volume. Plotted are the data from Fig. |S] (short-loop 
and pressure sensor), in comparison data using a smaller inter- 
droplet distance (pressure sensor and large loop, A ~ 3 mm). 
The empty symbol for the pressure sensor corresponds to a 
distance of A = 2.6 mm. The other parameters are rjd = 
1.4Pas, 77$ =0.02 Pas, Utot =6mm/s. 



dependence of the parameters is found when varying two 
of them at the same time. Figure [TT] shows such an ex- 
ample, namely the resistance length L^ as a function 
of the drop volume. The small-loop data have already 
presented in Fig. [Sj The corresponding large-loop resis- 
tances are systematically smaller. The error bars and an 
independent measurement with the pressure sensor both 
indicate that the difference is not an artifact. The dis- 
crepancy is misleading, since the difference here is due to 
the smaller inter-droplet distance of approximately 3 mm 
in the large-loop experiment. For the volume of 65 nl, we 
find the same value w 2 mm in Fig. [TO] at this very 
distance. The resistance length, when measured with the 
small-loop device, is approximately 5 mm at a volume of 
65 nl. This value is approximately the one of the plateau 
at large distances in Fig. [TOl We may thus conclude that 
the small-loop device really measures the resistances of 
individual droplets, and that this resistance corresponds 
to the plateau at large distances in the large-loop device. 

Figure[TT]leads us to a comment which we consider cen- 
tral in the analysis of droplet traffic in networks. When 
droplets are driven through a complex and possibly large 
network of microfluidic channels, comprising a lot of junc- 
tions, then no direct control over parameters such as their 
distance or their velocity is possible. When the droplet 
trains in an existing device exhibit irregularities, i.e. that 
they do not follow the trajectory they are supposed to fol- 
low in order to fulfill the objective of the device, then it 
will become necessary to understand the droplet dynam- 
ics in more detail. The present work demonstrates some 
aspects of this understanding for a simple device, where 
the mutual infiuence of the hydrodynamic resistance and 
the droplet movement becomes apparent. 
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V. SUMMARY 

In the present paper we use three methods to measure 
the effective hydrodynamic resistance of droplets in mi- 
crochannels: One standard method based on the use of 
pressure sensors and two variants of a new method re- 
lying on the analysis of droplet trajectories in a simple 
microfluidic device. These two variants are proposed as 
complementary techniques, as they are based on differ- 
ent model assumptions on the droplet behavior in the 
microfluidic network. 

The main idea put forward here is to use the 
global droplet dynamics, namely the distances between 
droplets, their number in each branch of the loop, and the 
choices taken at the first junction. All of them can simply 
be extracted by analyzing a video of the droplets flowing. 
The methods are non-invasive as they do not require ad- 
ditional measuring channels for pressure sensors, nor do 
they require expensive equipment. The interpretation of 
these observables to obtain the resistances is based on a 
theoretical model, the simplifying assumptions of which 
we test in the present work. The proposed methods focus 
rather on the understanding of the droplet dynamics than 
on the isolated measurement of the all-parameter depen- 
dence of the hydrodynamic resistance of droplets. In this 
sense, it is not meant to replace the standard method em- 
ploying pressure sensors, but rather to provide a guided 
example for how the droplet dynamics reacts on varia- 
tions of the hydrodynamic resistance and vice versa. 

We were able to show that there exists a parameter 
regime where the assumptions adopted by the model do 
hold. In this regime, the quantity that we extracted 
from the experimental observations indeed serves as the 
wanted hydrodynamic resistance of droplets. These re- 
sults are presented in Figs. O [HI El UHl and [TTl where 
we varied different parameters, namely the volume, ve- 
locity and viscosity of droplets, as well as their distance. 



All these sets of measurements were done in a way such 
that they avoid those parameter regimes where the model 
assumptions are not valid. Therefore, they do not con- 
tradict each other. This claim is enforced by independent 
and more direct measures of the droplet resistances, un- 
dertaken with an expensive and elaborate pressure-sensor 
technique. The requirement to find the (formerly un- 
known) valid part of a high-dimensional parameter space 
presented the main intellectual difficulty while generat- 
ing the results displayed in the figures. The parameter 
dependencies we found are the following: The hydrody- 
namic resistance increases with increasing volume and 
increasing viscosity of the droplet. It increases with inter- 
droplet distances as long as the distances are smaller than 
approx. 6 mm. Beyond that value the resistance does not 
depend on the inter-droplet distance. For small droplet 
velocities, the resistance strongly decreases with increas- 
ing velocity. It changes only little for velocities beyond 
3mm/s. On a more fundamental level of description, 
the results presented here demand a systematic deter- 
mination of the hydrodynamic resistance as a function 
of droplet volume, spacing, velocity and viscosity. This 
could be done for example by extending Refs. [ll, [2^ . 
which are concerned with the resistance of inviscid bub- 
bles, not droplets. In the limit of high velocities, large 
distances, and small volumes, the resistance will still have 
to turn out to be constant, as assumed by the simple 
model used here. 
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